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Transport and non-equilibrium magnetization in monolayers of magnetic molecules subject to a 
bias voltage are considered. We apply a master-equation approach going beyond the sequential- 
tunneling approximation to study the Coulomb-blockade regime. While the current is very small in 
this case, the magnetization shows changes of the order of the saturation magnetization for small 
variations of the bias voltage. Inelastic cotunneling processes manifest themselves as differential- 
conductance steps, which are accompanied by much larger changes in the magnetization. In ad- 
dition, the magnetization in the Coulomb-blockade regime exhibits strong signatures of sequential 
tunneUng processes c/e-exciting molecular states populated by inelastic cotunneUng. We also consider 
the case of a single molecule, finding that cotunneling processes lead to the occurrence of magnetic 
sidebands below the Coulomb-blockade threshold. In the context of molecular electronics, we study 
how additional spin relaxation suppresses the fine structure in transport and magnetization. 

PACS numbers: 73.63.-b, 75.50.Xx, 85.65.+h, 73.23.Hk 



I. INTRODUCTION 

The idea of molecular spintronics consists of integrat- 
ing the promising concepts of molecular electronics and 
spintronics ii^i^i^i^i^ A particularly interesting aspect of 
molecular electronics, besides the prospect of further 
miniaturization, is the possibility to use chemical syn- 
thesis for the fabrication of device components. This 
bottom-up process would start from relatively simply 
molecules and be massively parallel. In this context, 
spintronics is discussed in relation to magnetic memory^ 
and quantum computation*^ Both ideas rely on magnetic 
molecules Partly for this reason electronic transport 
through magnetic molecules has recently received a lot 
of attention . -^^•-^-^•-^^•-^^'-^'^•-^^•-^^•-^'^•-^^•-^^•^^•^-^•^^•^^•^^ 

Experimental research has focused on the fine struc- 
ture of the Coulomb-blockade peak^iiiii^iii and on 
Kondo correlations in single-molecule transistors ^iSiliiSi 
Furthermore, novel spin-blockade mechanisms and neg- 
ative differential conductance have been observed J^a^^ 
These findings have also stimulated theoretical work, 
which mostly employs the sequential-tunneling 
approximationJ^ii^ii^i2aiSii2^i2ii2^ Like artificial 
quantum dots, a molecular junction is in the Coulomb- 
blockade regime at sufficiently small bias voltage, 
except at crossing points where two states with electron 
numbers differing by unity become degenerate. Due 
to the discreteness of molecular many-particle energies 
for weak coupling to the leads, there are typically no 
molecular transition energies in the window between the 
chemical potentials of the leads for small bias. In this 
regime the very small tunneling current is not correctly 
described by the sequential-tunneling approximation. 
It is instead dominated by cotunneling^ which appears 
in fourth order in the perturbation expansion in the 
tunneling amplitude. However, despite its experimental 
observation^ transport through magnetic molecules in 
this regime has been little studied*^ 



We study magnetic molecules under a bias voltage in 
the Coulomb-blockade regime. Our main result is that 
while any features in the differential conductance are 
very small due to the suppression of the current, there 
are large changes in the average magnetic moments of 
the molecules with bias voltage and applied field. The 
measurement of magnetic moments of sub-monolayers of 
molecules has been demonstrated 20 years ago^ Even 
the detection of the spin of a single molecule may be 
feasible j^i^ However, it is not clear how to perform such 
a measurement in a molecular-junction experiment. One 
recent experiment suggests that it is possible to employ 
carbon nanotube superconducting quantum interference 
devices for the detection of the switching of singe mag- 
netic moments*^ We here mainly consider a monolayer 
of magnetic molecules between metallic electrodes, since 
the measurement of the magnetization of a thin film is 
expected to be easier than of a single molecule. Vari- 
ous molecules form nearly perfect monolayers on metallic 
substrates*^ 



To find the current and the non-equilibrium magnetiza- 
tion, we use the master-equation formalism, treating the 
tunneling to the leads as a perturbation ] ^^i^^i^^i^^i^^i^^ 
This approach describes the Coulomb and exchange in- 
teractions on the molecule exactly and works also far from 
equilibrium. In particular, it is not restricted to the li- 
near-response regime of small bias voltage. 



For memory applications the control of spin relaxation 
is crucial. Since cotunneling and additional spin relax- 
ation due to, e.g., dipolar and hyperfine interactions, 
have similar selection rules for molecular transitions, con- 
sistency requires to include both. We find that spin re- 
laxation is very effective in washing out the fine structure 
in the Coulomb-blockade regime but may be used to ad- 
vantage for the generation of spin-polarized currents. 
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FIG. 1: (Color online) Sketch of the geometry. A monolayer 
of magnetic molecules is adsorbed on a metallic substrate, 
which serves as a bottom electrode. A thin metallic layer is 
used as a top electrode. 



II. MODEL AND METHODS 



For the most part, we consider a monolayer of mag- 
netic molecules sandwiched between two metallic elec- 
trodes, see Fig. [H We assume that magnetic interac- 
tions between the molecules are negligible and that all 
molecules have the same spatial orientation relative to 
the electrodes. In this case it is sufficient to consider 
the properties of a single molecule. Relaxation in the 
leads is assumed to be fast so that their electron distribu- 
tions can be described by equilibrium Fermi functions. In 
the simplest case, transport involves tunneling through 
only a single molecular level with onsite energy Cd and 
local Coulomb repulsion U. The full Hamiltonian of the 
system reads H = i^moi + ^leads + Ht , where^ 



U 



— rid {rid - 1) - Js • S 



(1) 



describes the molecular degrees of freedom, i^ieads = 
Z)a=L,R Zlka eakaj,ka^«ka represents the two leads a = 
L,R (left, right), and Ht = Ea=L,R ^« Eka(«LaC^ + 
cj- acker) describes the tunneling. The tunneling ampli- 
tudes ta are chosen real. The operator cj. creates an 



electron with spin a on the molecule, rid ^ 



and s = Xl^^/ c'laaa'Ca' /2 are the corresponding num- 
ber and spin operator, respectively. The parameter J 
denotes the exchange interaction between the electrons 
and a local spin S, where S ■ S = S{S -\- 1). We restrict 
ourselves to the case of easy- axis anisotropy, K2 > 0. For 
simplicity we consider identical g factors for s and S. An 
external magnetic field B is applied along the easy axis 
of the molecule, where a factor gjiiB has been absorbed 
into B. aj^j^^ creates an electron in lead a with spin a, 
momentum k and energy e^h- 

The leading contribution to the transition rates be- 
tween molecular many-particle states is of second order 
in i^t, corresponding to sequential tunneling. The tran- 



sition rates can be obtained from Fermi's Golden Rule^^ 

a 

+ [l-f{e^-e^-^ia)]\CU^). (2) 

Here, the eigenstates \n) and \n') of i^moi denote the 
initial and final state of the molecule, respectively, u^a 
is the density of states per unit cell of electrons with 
spin a in lead a, /(e) is the Fermi function, fia is the 
chemical potential in lead a, where /iL — Mr = — eF, and 
C^^, = {n\C(j\n') is the matrix element of the electron 
annihilation operator between molecular many-particle 
states. The typical sequential-tunneling rate involving 
lead a and electrons with spin a is given by l/r^a = 
27r t^^Uaa is sct to unity). 

To go beyond the leading order, the tunneling Hamil- 
tonian is replaced by the T matrix^ which is self- 
consistently given by 



Ht + Ht 



1 



Ei - Ho^if] 



T. 



(3) 



Here, Ei is the energy of the initial state |i)|n), where \i) 
refers to the equilibrium state of the left and right leads 
(at different chemical potential) and \n) is a molecular 
state. Furthermore, Hq = i^moi + ^ieads, with the energy 
of the leads measured relative to equilibrium, and 7^ is a 
positive infinitesimal ensuring that the Green function in 
T is retarded. To fourth order, the transition rate from 
state 1^) |n) to |/) \n') with an electron tunneling from lead 
a to lead a' is given by 
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The energies of the initial state \n)\i) and final state 

1^01/) = I^O^Lk'a'^akaK) are denoted by Ei and Ef, 
respectively. We restrict ourselves to the case of infinite 
/7, i.e., double occupancy of the molecule is forbidden. 
Inserting Ht and summing over final lead states yields 
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where T^^/^^ ^^''^^) denotes the cotunneling rate de- 
scribing virtual transitions between two empty (singly oc- 
cupied) molecular states. Here, we have assumed that the 



density of states in the leads is independent of energy. To 
the same order in i^t, one also obtains processes chang- 
ing the electron number by ±2. For U ^ oo these pair- 
tunneling processes^^ are suppressed. Note that Eqs. ([5|) 
and dn]) contain both elastic and inelastic cotunneling. 

Since the above expressions diverge due to second- 
order poles from the energy denominators, the cotun- 
neling rates cannot be evaluated directly i^^i^^i^^i^^ We 
apply a regular izat ion scheme that follows Refs. [38||39||4Q| 
and is motivated by the observation that Eqs. ([5]) and 
(|6]) do not take into account that the intermediate state 
obtains a finite width F due to tunneling. In our regime 
of weak tunneling, the width F is of second order in the 
tunneling amplitudes ta- This width is introduced into 
the energy denominators, replacing r]. When the cotun- 
neling rates are expanded in powers of F, it turns out 
that the leading term is of order 1/F ex This can- 

cels two powers of the tunneling amplitude in Eqs. ([5]) 
and (|6]) so that the result is in fact a sequential-tunneling 
contribution. Since we have already included the full 
sequential- tunneling rates, this new contribution should 
be dropped. We thus take the next order, F^, for the 
cotunneling rates. 

The ad-hoc regularization of the cotunneling rate is not 
necessary in a description of cotunneling through wide 
quantum dots using non-equilibrium Green functions.— 
This approach avoids the divergences and also leads to 
a renormalization of transition energies at fourth order 
in the tunneling Hamiltonian. For a fully quantitative 
description of cotunneling through magnetic molecules it 
would be desirable to employ this approach, which is, 
however, made complicated by the presence of the in- 
ternal spin degree of freedom. On the other hand, the 
T-matrix approach used here is found to give qualita- 
tively reasonable results in comparison with cotunnel- 
ing experiment^ and we expect it to catch the relevant 
physics for the system studied here. 

The sequential and cotunneling rates appear in the rate 
equations for the probabilities to find the molecule in 
state |n) (we assume rapid dephasingi^) , 
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non-zero only if both \n) and \m) are empty (singly oc- 
cupied). The current through the left lead is given by 
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The steady-state probabilities of molecular states are 
obtained by solving Eq. ([7j) with the time derivatives set 
to zero. The average magnetization in the z direction per 
molecule is given by M = rrinP^, where rrin denotes 
the quantum number of the z component of the total spin 
s + S in state |n). 




3eJ4e 

d 



FIG. 2: (Golor online) (a) Differential conductance dl/dV 
and (b) probabilities of molecular many-particle states 
as functions of bias voltage V, for low bias voltages. The 
probability P^/^ of the ground state has been scaled by a 
factor 1/3. Here, we assume S — 2^ J — K2 — 5meV, td — 
10 J, B = 2 meV, and T = 0.3 meV.i^ 



III. RESULTS AND DISCUSSION 

We start by discussing the results obtained for the 
differential conductance dl/dV at low bias voltages. If 
the system is in the Coulomb-blockade regime, sequen- 
tial tunneling is thermally suppressed and transport is 
dominated by cotunneling. The magnitude of the cur- 
rent is then small. The conductance at zero bias voltage 
is finite, see Fig. [2](a), due to elastic cotunneling. The 
cotunneling rates are proportional to the bias voltage, if 
the molecular level is far from the chemical potentials, 
leading to ohmic behavior. The rounded steps in dl/dV 
correspond to the onset of additional inelastic cotunnel- 
ing processes. Selection rules for the spin quantum num- 
ber require Am = 0, ±1. For the parameters chosen in 
Fig. [21 the ground state has electron number n = 1 and 
maximum spin, m = 5/2. Inelastic cotunneling processes 
corresponding to the two steps involve the two different 
final states with n = l, m = 3/2 and virtual occupation 
of the state with n = 0, m = 2, as illustrated in Fig. [3l 
Further steps in dl/dV are not observed, since the cor- 
responding inelastic cotunneling transitions have smaller 
energy differences between initial and final states and are 
therefore activated immediately when the probability of 
the initial state becomes significant. 

Cotunneling steps and sequential tunneling peaks show 
fundamentally different dependences on the onsite energy 
Cd. For single-molecule junctions it is possible to change 
Cd by applying a gate voltage, e.g., in molecular-junction 
experiments. For monolayers one does not have this op- 
portunity. We come back to this point below. While the 
bias voltages at which sequential tunneling peaks occur 
shift linearly with e^^, the positions of cotunneling steps 
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FIG. 3: (Color online) Level scheme showing the energies of 
molecular states as a function of magnetic quantum number 
m for electron numbers n = 0, 1.— The dashed double arrows 
signify inelastic cotunneling between the ground state with 
m — 5/2 and the two states with m = 3/2, involving vir- 
tual occupation of the state with n = and m = 2. While 
sequential tunneling requires a change of the electron num- 
ber by An = dzl and of the magnetic quantum number by 
Am = ±1/2, cotunneUng processes obey the selection rules 
An = 0, Am = 0,±1. 
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remain unaffected. This follows directly from evaluat- 
ing Eqs. ©-(in]) in the limit of large Cd-^^ For magnetic 
molecules, the position of the cotunneling steps shifts lin- 
early as a function of the external magnetic field due to 
the Zeeman effect, as observed for Mni2.^^ 

While dl/dV represents the change of the very small 
current with bias voltage in the cotunneling regime, the 
change of the probabilities of molecular states with 
bias voltage is of order unity, as shown in Fig.[2](b). The 
probability of the lowest-energy state with m = 5/2 de- 
creases, whereas the probabilities of other states increase. 
Cotunneling enables transitions between molecular states 
with the same electron number but with magnetic quan- 
tum numbers differing by Am = ±1. These transitions 
are suppressed only as the inverse square of the energy 
difference between the initial state and the virtual state 
involved. In sequential tunneling, such transitions are 
also possible, requiring two consecutive steps, but are ex- 
ponentially suppressed in the Coulomb-blockade regime. 
In the sequential-tunneling approximation the molecule 
would thus remain in the lowest-energy state with essen- 
tially unit probability. This approximation is evidently 
invalid for determining the probabilities in this regime. 

Interestingly, the strong effect of cotunneling on the 
probabilities also leads to observable effects of sequential 
tunneling on transport in the cotunneling regime. ^^^^^ 
While sequential tunneling starting from the lowest- 
energy state is exponentially suppressed, sequential tun- 
neling from higher-energy states can be possible. With 
increasing bias voltage, these higher-energy states be- 
come increasingly populated due to cotunneling, as 
Fig. WO^) shows. This leads to sidebands in dl/dV in 
the Coulomb-blockade regime that show the linear de- 
pendence on the gate voltage characteristic of sequen- 
tial tunneling^ Strong electron-phonon coupling can en- 



FIG. 4: (Color online) (a) Magnetization M, (b) linear plot of 
the current /, (c) logarithmic plot of the sequential-tunneling 
current, and (d) probabilities of various molecular many- 
particle states as functions of bias voltage V. The parameters 
are chosen as in Fig. [2] 



hance this effect, since it crucially affects the ratio of the 
rates for sequential and cotunneling processes^Si In our 
case, these sidebands are very weak, since the current is 
controlled by the small cotunneling rates. However, we 
will see that the effect on the probabilities P^ of molec- 
ular states is significant. 

Figure [H^a) shows the average magnetization per mo- 
lecule as a function of bias voltage over a broad range 
including both the cotunneling and sequential-tunneling 
regimes. The magnetization is nonzero due to an external 
magnetic field. At zero bias, the molecule is in its ground 
state with m = 5/2. The onset of inelastic cotunneling 
to the two states with m = 3/2 leads to a decrease in the 
magnetization in each case. 

The bias-voltage dependence of the magnetization for 
voltages above the Coulomb blockade threshold is ac- 
companied by sizeable steps in the current, as seen in 
Fig. mb). At each of these fine-structure steps an addi- 
tional inelastic sequential-tunneling transition becomes 
possible. The Coulomb-blockade threshold corresponds 
to the transition with initial state n = l, m = 5/2 and 
final state n = 0, m = 2. Therefore, the onset of se- 
quential tunneling is accompanied by a decrease in the 
magnetization. At large bias the magnetization drops to 
zero since all states are occupied with equal probability. 

Remarkably, pronounced step-like features are also 
present below the Coulomb-blockade threshold in 
Fig. SJa), where the current is due to cotunneling and 
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FIG. 5: (Color online) Level scheme illustrating the interplay 
between sequential tunneling (solid arrows) and cotunneling 
(dashed arrows) in magnetic molecules.— Even below the 
Coulomb-blockade threshold sequential tunneling processes 
starting from higher-energy states populated by cotunneling 
may cause the depopulation of these states and drastically 
affect the average magnetization. At the step denoted by 
an arrow in Fig. HJa), the excitation of the transition with 
m=— 1/2^0 (heavy solid arrow) gives rise to a redistribu- 
tion of the probabilities . Note that exothermal transitions 
with Am = ±1/2 are always possible. 



thus very small, cf. Figs. [H^b), (c). This can be under- 
stood from the bias-voltage dependence of the relevant 
probabilities in Fig. W^d). As an example, consider 
the step marked by an arrow in Fig. [H^a). The physics 
leading to the drastic change of the probabilities is illus- 
trated in Fig. [5l The sequential tunneling processes with 
m = -3/2 ^ -2, m = -1/2 ^ -1, m = 3/2 ^ 2, 
and m = 1/2 ^ 1, starting at the higher-energy level of 
each pair (thin arrows in Fig. [5j), are already energeti- 
cally possible at lower bias voltages causing the partial 
depopulation of the initial states. However, the prob- 
abilities of these states are non-zero mainly due to co- 
tunneling processes (dashed arrows in Fig. [5]). Below the 
step marked in Fig. [41(a), the half-integer spin states with 
positive and negative m are not connected by sequen- 
tial tunneling processes. As soon as the transition with 
m = — 1/2^0 (bold arrow in Fig. [5]) becomes possible, 
the states with positive and negative m are connected and 
fast, sequential tunneling processes depopulate all states 
except for the ground state, which has m = 5/2. Con- 
sequently, the average magnetization again approaches 
its maximum value. Similarly, one can attribute each 
step to a particular molecular transition. As Fig. HJc) 
shows, the onsets of some of these sequential-tunneling 
processes can also be seen in the sequential-tunneling cur- 
rent, which is, however, tiny in the cotunneling regime. 

The above discussion shows that quantities that de- 
pend strongly on the probabilities of molecular states, 
such as the magnetization, are much more sensitive to 
changes of the bias voltage in the Coulomb-blockade 
regime than the conductivity. This suggests to use 
the magnetization-voltage characteristics, i.e., the mag- 
netization as a function of bias voltage, instead of the 
current-voltage characteristics to extract the excitation 
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FIG. 6: (Color online) Magnetization M as a function of bias 
voltage V and magnetic field B for different spin relaxation 
times: (a) trei = oo, (b) trei = lO^r, (c) trei = lO^T, and (d) 
trei — 0. Here r = {27Ttai^a)~'^ denotes the typical electronic 
tunneling time, assuming symmetric coupling to the leads. 
All other parameters are chosen as above. The dashed lines 
denote the Coulomb-blockade threshold, (e) Level schemes 
illustrating the origin of the magnetization plateaus. 



spectrum of magnetic molecules. In order to distinguish 
magnetic transitions from, e.g., vibrational excitations, 
one should analyze their dependence on the magnetic 
field. Furthermore, for a monolayer there is no gate volt- 
age that can serve as an independent parameter. The 
magnetic field can assume this role. 

Figure [6](a) shows a density plot of the magnetization 
as a function of bias voltage and magnetic field. The mag- 
netization is an odd function of the field. The transition 
energies shift linearly with the field, = Am 5, if the 
initial and final states have magnetic quantum numbers 
differing by Am. 

Complementary to conventional differential-conduc- 
tance plots, the density plots in Fig. [6] can serve as finger- 
prints of the internal degrees of freedom of the molecules. 
The Zeeman splitting of the molecular levels due to the 
external magnetic field gives rise to triangular plateaus 
with a tip at 5 = 0. These plateaus are bounded by two 
sequential- tunneling transitions. In each case, these two 
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transitions differ in the sign of the magnetic quantum 
number m of both initial and final molecular states. For 
the chosen parameters, the plateaus can be attributed 
to the following transitions from empty to singly oc- 
cupied states, starting at low bias voltage (cf. Fig. [5]): 
\m\ = 3/2 ^ 2, |m| = 1/2 ^ 1, |m| = 1/2 ^ 0, 
\m\ = 3/2 ^ 2, |m| = 3/2 ^ 1, \m\ = 1/2 ^ 1, 
\m\ = 5/2 2 (this is the first transition starting 
from the ground state and thus represents the Coulomb- 
blockade threshold), |m| = 1/2^0, and \m\ = 3/2 1. 
Several transitions appear twice because there are two 
states with magnetic quantum numbers ±3/2 and ±1/2, 
respectively. For a local spin S = 2 there exist nine tran- 
sitions obeying the selection rule Am = ±1/2, as can be 
seen from Fig. [5l in accordance with the nine plateaus 
shown in Fig. [61(a). Note again that the signal is similar 
on both sides of the Coulomb-blockade threshold. 

The origin of the plateaus is schematically illustrated 
in Fig. [61(e) for the transition \m\ = 1/2 0. In the ab- 
sence of an external Zeeman field the excitation energies 
for both transitions is equal. However, the excitation en- 
ergies differ by the Zeeman energy as soon as a magnetic 
field is switched on. This leads to the occurrence of a 
finite bias- voltage window, where the excitation of one of 
the two transitions is energetically possible whereas the 
other one is not. Inside this window, only the spin-down 
state is depopulated by sequential tunneling, leading to 
a large positive magnetization. 

So far we have restricted ourselves to the situation 
where the relaxation of the local molecular spin is dom- 
inated by electron tunneling, i.e., the spin is conserved 
between tunneling events. However, there are other pro- 
cesses that also contribute to spin relaxation: (i) Mag- 
netic molecules containing transition-metal ions, such as 
Mni2 clusters, show strong spin-orbit interaction, which 
leads to spin relaxation, (ii) Hyperfine interactions with 
nuclear magnetic moments in the molecule can also lead 
to spin relaxation. However, in molecules one has the 
chance to essentially remove this mechanism by choosing 
isotopes with vanishing nuclear spins, (iii) Dipolar inter- 
actions with spins of other molecules in the monolayer or 
with impurity spins in the electrodes contribute to spin 
relaxation, (iv) Small non-uniaxial magnetic anisotropics 
lead to tunneling between the eigenstates of i^moi- This 
mechanism has recently been discussed in the context of 
transport through magnetic molecules J^ii^ii^ 

All these processes change the magnetic quantum num- 
ber while keeping the electron number constant (An = 
0). The dominant transitions are the ones with Am = 
±1. These are the same selection rules as for cotunnel- 
ing, indicating that one should include additional spin 
relaxation for consistency when studying cotunneling. 

The effect of spin relaxation on the electronic transport 
is included in the formalism by a phenomenological rate 
(X 1 / trei which forces the system to approach the equilib- 
rium distribution on the timescale trei- We include ad- 
ditional transition rates between states \n) and |m) with 
An = and Am = ±1, F^^^ = exp[(en - 6^)//cT]/trei 
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FIG. 7: (Color online) (a), (b) Magnetization M and (c), (d) 
differential conductance of a single magnetic molecule as a 
function of V and ea for (a), (c) slow spin relaxation, trei = 
lO^V, and (b), (d) fast spin relaxation, trei = t. We assume 
5' = 2, J = ^2 = 5meV, T = 0.1 meV, and B = 2meV. 



for Cn < Cm and FJ^f^ = l/ty-ei otherwise. The additional 
rates obey detailed balance, ensuring relaxation towards 
equilibrium in the absence of tunneling. 

Effects of spin relaxation on the bias-voltage depen- 
dence of the magnetization are illustrated in Figs. [H^a)- 
(d). For small trei (fast relaxation), the number of tran- 
sitions appearing as steps in the magnetization-voltage 
characteristics is reduced, since spin relaxation depopu- 
lates higher-energy states that serve as initial states for 
these transitions. 

The magnetization plateaus start to occur when the 
relaxation time trei becomes significantly larger than 
the typical sequential-tunneling time r = (27rt^z^c^)~^. 
(The sequential-tunneling time enters because the rele- 
vant process is the depopulation of states by sequential 
tunneling.) Then the time spent by the electron on the 
molecule is smaller than the spin relaxation time so that 
magnetic excitations survive between tunneling events. 

So far we have considered a monolayer of magnetic 
molecules, mostly because the measurement of the mag- 
netization is easier for larger numbers of molecules. 
As mentioned previously, even the detection of a sin- 
gle molecular spin might be feasible Using a single 
molecule allows one to introduce a gate electrode in or- 
der to tune the molecular energy levels by shifting e^^, see 
Eq. ([1]). In the following, we briefly discuss results ob- 
tained for varying gate voltage. To increase the magneti- 
zation signal while retaining the gate electrode, one might 
consider a one-dimensional array of magnetic molecules 
or even a large number of such arrays aligned in parallel. 

The plot of the magnetization and the differential con- 
ductance as functions of bias voltage and onsite energy 
Cd presented in Figs.[7Ka),(c) shows two striking features. 
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FIG. 8: (Color online) Polarization of the current, p = (/^^ — 
/^^)/(/^^ + /^^), as a function of spin relaxation time trei in 
units of the typical tunneling time r. 



First, the magnetization shows steps indicating the onset 
of inelastic cotunneling which are almost independent of 
€d. The corresponding steps in dl/dV are very small in 
absolute units, see Fig. [2](a). 

Second, the magnetization shows strong additional 
magnetic sidebands in the Coulomb-blockade regime. 
These sidebands are the consequence of sequential- 
tunneling transitions depopulating molecular states that 
are populated by cotunneling, as discussed above. In 
dl/dV the corresponding features are completely hidden 
by the low-bias tail of the large peak at the Coulomb- 
blockade threshold (not shown). The observation of these 
sidebands in the Coulomb-blockade regime requires spin 
relaxation times long compared to the typical tunnel- 
ing time. For fast spin relaxation, fine-structure peaks 
are only present in the sequential- tunneling regime, see 
Fig. [71(b), since sequential tunneling is still faster than 
spin relaxation, even though cotunneling is slower. As 
shown in Fig. E^d), the absence of such sidebands is 
accompanied by suppressed fine-structure peaks in the 
sequential-tunneling regime. 

Finally, we note that sufficiently fast spin relaxation 
leads to spin-polarized stationary currents in the pres- 
ence of a magnetic field. If the spin of the magnetic 
molecule relaxes fast compared to the typical tunneling 
rate, which is essentially determined by the current, the 
system is always in its ground state. Due to the Zeeman 
effect the ground state has maximum magnetic quantum 
number; m = 5/2 for our example. Thus only spin-down 
electrons can tunnel onto the molecule, resulting in a 
spin-polarized current. Note that this argument is not 
restricted to low-order perturbation theory in Ht. As 
shown in Fig. [51 the degree of spin polarization is basi- 
cally determined by the ratio of the spin relaxation rate 



and the typical electronic tunneling rate. 

IV. CONCLUSIONS 

In summary, we have studied the interplay of electronic 
transport through magnetic molecules and their non- 
equilibrium magnetic moment beyond the sequential- 
tunneling approximation. We have focused mostly on 
monolayers, which should give a better chance to mea- 
sure the magnetization than single molecules would. 

While the excitation of inelastic tunneling processes in 
the Coulomb-blockade regime leads only to a very small 
absolute change in the current, the change of the prob- 
abilities to find the molecule in various many-particle 
states is significant. This manifests itself in a strong bias- 
voltage dependence of the magnetization. The magneti- 
zation of a molecular monolayer can be switched by an 
amount of the order of the saturation magnetization by 
a small change of bias voltage, and without causing the 
flow of a large current. 

We find steps in the differential conductance due to 
inelastic cotunneling, which have been observed in ex- 
periments on Mni2.^^ These steps are accompanied by 
much larger changes in the magnetization. Another in- 
teresting effect is the appearance of additional sidebands 
in the Coulomb-blockade regime that can be ascribed to 
(ie-excitations by sequential tunneling of states populated 
by cotunneling. These sidebands are very prominent in 
the magnetization. We suggest that the magnetization, 
or any measurable quantity that strongly differs between 
molecular states, can be employed to study molecular 
transitions that are, from the point of view of transport, 
hidden in the Coulomb-blockade regime. 

For spintronics applications, the ability to control 
the persistence of the stored information is crucial. In 
this context, we have considered effects of additional 
spin relaxation in the same formalism. Our results 
show that for sufficiently fast spin relaxation the peaks 
in the differential conductance and the steps in the 
magnetization are washed out, as expected. At the 
same time, the degree of polarization of the steady-state 
current contains information about the ratio of the spin 
relaxation rate and the typical electronic tunneling rate. 
Fast spin relaxation, while in general undesirable, can 
lead to a highly polarized current in the presence of a 
magnetic field. 
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